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摘 要 : 本 文 围绕 统计 过 程 的 监测 、 诊 断 和 变 点 时 间 估 计 三 个 基本 问题 ， 首 先 概括 地 介绍 了 过 程 监测 
的 最 优 控制 图 一 累积 和 控制 图 ; 然后 ， 详 细 介 绍 了 如 何 利用 多 重 累积 和 控制 图 对 过 程 的 异 变 
进行 统计 诊断 ;最 后 简 述 了 三 种 常用 的 变 点 估计 方法 : 最 大 似 然 估 计 ，Schwarz 信息 判别 估计 
和 Bayesian 估计 。 此 外 ， 文 中 还 列 出 6 个 囊 待 解决 的 问题 。 

关键 词 : 变 点 监测 ， 统 计 诊 断 ， 质 量 管理 ， 统 计 过 程控 制 

分 类 号 : AMS(2000) 62L10 中 图 分 类 号 : 0213.1; 0212.3 文献 标识 码 : A 


1 引言 


我 们 知道 ， 无 论 一 个 产品 的 生产 过 程 ， 还 是 一 个 服务 系统 的 运行 ， 都 无 不 受到 各 种 因素 的 
干扰 。 统 计 过 程控 制 (Statistical Process Control, SPC) 就 是 应 用 统计 分 析 技 术 对 所 考察 的 系 
统 运行 过 程 进行 监控 ， 科 学 地 区 分 出 影响 过 程 正 常 运行 的 随机 因素 (不 可 避免 的 因素 ) 和 异常 因 
素 ， 进 而 对 过 程 发 生 异 常 变化 的 时 间 、 大 小 、 范 围 ， 甚 至 原因 提出 警示 ， 以 便 及 时 采取 措施 ， 
消除 异常 ， 使 过 程 处 于 正常 的 运行 状态 。 

SPC 的 基本 原理 和 方法 是 上 世纪 30 年代 由 Shewhart 博 士 为 了 有 效 地 对 生产 过 程 中 产品 质 
量 进行 监测 控制 而 提出 的 ， 至 今 已 有 70 多 年 的 历史 。 自 创立 以 来 ， 它 就 在 工业 和 服务 等 行业 
得 到 了 推广 和 使 用 。 二 战 时 期 美国 将 其 制定 为 战 时 质量 管理 标准 ， 为 保证 军工 产品 的 质量 和 及 
时 交付 起 到 了 重要 作用 。 战 后 的 日 本 从 1950-1980 年 在 工业 界 广泛 推广 和 应 用 SPC， 使 日 本 跃 
居 世 界 产品 质量 和 生产 率 的 领先 地 位 ， 以 至 于 美国 著名 的 质量 管理 专家 Berger 教授 也 曾 说 ; 日 
本 成 功 的 基石 之 一 就 是 SPCI。 从 上 世纪 80 年 代 起 ，SPC 在 许多 工业 发 达 国 家 复兴 ， 世 界 很 
多 大 公司 也 纷纷 在 自己 内 部 积极 推广 和 应 用 SPC。 虽 然 ， SPC 是 从 产品 的 质量 监控 开始 的 ， 但 
经 过 70 多 年 实践 和 发 展 ， 尤 其 是 与 计算 机 技术 的 紧密 结合 ， 其 原理 和 方法 现 已 广泛 应 用 于 设 
计 、 销 售 、 服 务 、 管 理 等 过 程 。 

SPC 的 广泛 应 用 大 大 促进 了 SPC 理论 研究 的 发 展 。 本 文 将 主要 围绕 SPC 的 三 个 基本 问题 : 
监测 、 诊 断 和 变 点 时 间 估 计 问 题 ， 概 括 地 介绍 它们 的 基本 概念 、 方 法 以 及 与 我 们 的 研究 工作 有 
关 的 一 些 新 的 结果 ， 并 列 出 若干 重 待 解决 的 问题 。 


2 监测 问题 


所 谓 SPC 的 监测 问题 ， 狭 义 上 理解 就 是 如 何 利用 所 获得 的 被 监测 过 程 的 统计 数据 ， 构 造 控 
制图 ， 使 之 当 所 监测 的 过 程 有 蜡 常 变化 时 能 迅速 地 报警 。 由 于 被 监测 过 程 大 都 是 在 随机 因素 
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影响 的 环境 中 进行 ， 因 而 我 们 可 将 其 看 成 是 随机 过 程 。 既 然 是 随机 过 程 ， 当 对 其 进行 监测 报 
警 时 ， 就 可 能 存在 误 报 问 题 。 因 此 ， 监 测 问题 的 关键 就 变 成 如 何 构造 最 优 的 控制 图 ， 使 之 当 
所 监测 的 随机 过 程 有 异常 变化 时 能 在 很 小 的 误 报 率 前 提 下 ， 最 快 地 报警 。 用 数学 符号 表述 如 
Be WA, t= 1,2,…， 是 从 所 监测 的 随机 过 程 中 得 到 的 样本 (统计 数据 )， 其 中 Xe 表示 过 程 
在 t 时 刻 所 处 的 状态 。 设 在 时 刻 7 ( 称 为 变 点 时 刻 )， 过 程 发 生 了 异常 变化 ， 也 即 ， 过 程 在 + 时 
刻 之 前 服从 概率 分 布 pp， 而 在 7 时 刻 之 后 变 为 分 布 p1。 为 了 检测 这 一 变化 ， 我 们 可 以 利用 所 得 
到 的 统计 数据 ， 并 对 其 进行 加 工 和 处 理 ， 也 即 ， 选 取 适 当 的 统计 量 ，f.(X1,X2,.… ,Xn)， 其 
P fS fallt ,Zn) 是 变量 Z1,72,… ,zn 的 n 元 函数 。 进 而 ， 构 造 控制 图 


T(f)= min{n>1: fn(X1, X2,° n) >c}, (1) 


其 中 ec > 0 是 上 控制 界线 。 实 际 上 ， 一 个 控制 图 就 定义 了 一 个 首次 报警 (过程 有 异常 变化 ) 的 时 
刻 ， 在 随机 过 程 理论 中 称 这 个 时 刻 为 停 时 。 显 然 ， 选 取 不 同 的 统计 量 ， 就 可 定义 不 同 的 控制 
图 。 所 谓 最 优 控制 图 可 定义 如 下 。 设 有 统计 量 f* 和 按 式 (1) 定义 的 控制 图 TT(f*) 满足 : 对 任意 
BHRT) = ETS) = ETEN, A 


Ex, (T(f*)) < Ei(T()), (2) 


我 们 称 这 样 的 控制 图 T(f*) 为 最 优 控制 图 ， 其 中 Eo(-) 和 Ei() 分 别 表示 由 概率 分 布 po 和 pi 所 
确定 的 数学 期 望 。 也 就 是 说 ， 在 同样 的 平均 误 报时 间 的 条 件 下 ， 最 优 控制 图 了 (了 *) 在 所 有 控制 
图 中 最 先 报警 。 

人 我 们 先 来 回顾 一 下 已 提出 和 研究 的 一 些 控制 图 。 自 1931 年 Sh- 
ewhart!? a 入 们 先后 提出 了 各 种 各 样 的 控制 图 。 如 ， 多 变量 过 
Fes tA), BAA (cumulative sum, CUSUM) 控制 图 和 由， 指数 加 权 移 动 平均 (exponentially 
weighted moving average, EWMA) 控制 图 问 ，Shiryayev-Roberts 控制 图 [56,14， 贝 叶 斯 统计 控制 
图 外， 回归 模型 控制 图 加，Cuscore 控 制图 上 0， 残 差 过 程控 制图 上 14， 广义 似 然 比 控制 图 D3， 
轮廓 控制 图 13， 比 例 积 微分 控制 图 上 4， 广 义 EWMA 控 制图 415 ， 无 参数 Cuscore 控制 图 116] 等 
等 。 

显然 ， 在 众多 的 控制 图 中 ， 寻 找 最 优 控制 图 并 非 易 事 。 当 取 统 计量 


fr(®1,%2,°°° , Ln) = max Diels rr) |， 


1<j<n 
JEP L(x) Æ pi 关于 po HY Randon-Nikodym FAK, WF (1) 式 定义 的 控制 图 


Wl aa 2c}, (3) 


就 是 在 前 面 已 提 到 的 由 Page 于 1954 年 提出 的 CUSUM 控制 图 。 
19714, Lorden”! 首先 证 明了 当 监 测 样本 {Xi,k > 1} 是 独立 同 分 布 的 随机 过 程 时 ， 在 下 
述 Lorden 意义 下 CUSUM 控制 图 是 渐进 最 优 的 ， 若 Bo(Tp) = ET FLA, M 


Ei(Tp) < Fi(T) £ supesssup EN{(T—7+1)t|X,.. Xr], (4) 
7 之 1 


其 中 工 是 任意 的 一 个 控制 图 ， 而 (MM(-) 表 示 在 7 时 刻 之 后 过 程 服从 分 布 pi 所 对 应 的 数学 期 
望 。 
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这 里 需要 说 明 的 是 ， 上 述 关 于 最 优 控制 图 的 两 种 定义 (2) 和 (4) 式 是 不 等 价 的， 而 且 (4) 式 
要 比 (2) 式 更 难于 验证 。1986 和 1990 年 ，Moustakidesllal 和 Ritov09 先后 用 不 同 的 方法 证 明 
T% Eo(Tp) WA BRIN, CUSUM 控制 图 在 Lorden 意义 下 仍然 是 最 优 的 ， 即 ， 满 足 (4) 式 。 最 
近 ， 我 们 进一步 证 明了 当 监 测 样本 {Xk,k > 1} 是 平稳 的 线性 过 程 时 ，CUSUM 控制 图 依然 是 
ERRAL WA, RARS N: 

问题 1 对 于 一 般 的 随机 过 程 ，CUSUM 控制 图 是 否 是 最 优 的 ? 

在 上 述 最 优 控 制图 的 介绍 中 ， 我 们 有 一 个 很 强 的 假设 : 过程 发 生 异 常 变化 后 ， 其 服从 的 概 
率 分 布 是 已 知 的 。 实 际 上 ， 当 所 监测 的 过 程 发 生 异 变 后 ， 我 们 很 难 准确 地 知道 它 变化 后 服 
从 什么 概率 分 布 ， 尤 其 是 在 线 监测 时 更 是 如 此 。 这 样 ， 我 们 就 有 问题 2。 

问题 2 ”如 果 只 知道 所 监测 随机 过 程 在 正常 运行 时 (在 过 程 发 生 异 常 变化 之 前 ) 所 服从 的 概 
率 分 布 ， 我 们 能 否 找到 最 优 的 控制 图 ? 

当然 ， 关 于 随机 过 程 的 监测 还 有 许多 其 它 的 问题 ， 这 里 就 不 再 涉及 了 。 


3 诊断 问题 


这 里 所 说 的 诊断 是 指 统计 诊断 ， 它 是 利用 被 监测 过 程 的 统计 数据 ， 对 过 程 异常 变化 的 大 
小 ， 程 度 ， 范 围 以 及 变化 的 方向 和 趋势 等 及 时 地 做 出 统计 推断 。 在 过 去 的 几 十 年 ， 通 过 构造 各 
种 控制 图 来 对 所 监测 过 程 的 异常 变化 进行 诊断 已 有 很 多 研究 工作 ， 这 可 参阅 文献 [21-29]。 下 面 
主要 介绍 我 们 最 近 两 年 的 一 个 研究 结果 B99， 以 此 站 明 如 何 通 过 控制 图 来 进行 诊断 。 

设 所 监测 的 随机 过 程 Xe, t = 1,2,:…， 在 正常 情况 下 (在 异常 变化 之 前 ) 服从 概率 分 布 po,， 
其 中 po 既 可 表示 分 布 函数 也 可 表示 分 布 密度 函数 ， 而 ao 是 一 个 参数 。 假 设 我 们 根据 监测 
对 象 的 有 关 专 业 知 识 可 以 知道 ， 所 监测 的 随机 过 程 发 生 异 常 变 化 后 (在 变 点 时 刻 r 之 后 )， 其 
所 有 可 能 服从 的 概率 分 布 po 构成 一 个 集合 {po,a E 时 ， 其 中 列 是 一 个 不 含 ao 的 参数 集合 ， 
并 且 满 足 : 4aad e AHa 4a t, p £ pa'。 上 述 假设 意 指 ， 我 们 仅 知 道 过程 异 常 变 
化 后 其 所 有 可 能 服从 的 概率 分 布 为 {pa, a E 对 ， 但 具体 服从 哪 一 个 概率 分 布 并 不 知道 。 简 
记 T( = {po, Qa E A}。 通 常 集合 T(20) 含有 无 穷 多 个 概率 分 布 。 比 如 ， 某 过 程 在 正常 情况 下 
服从 标准 正 态 分 布 N(0,1)， 即 


一 1 2 
Poo = pol(z) = (V2r) ez /2, 


这 里 ao = (0,1)， 而 在 异常 变化 后 它 可 能 服从 的 分 布 为 


{Pa = Pino) (£) : p1 SMS pea, WF0;0< 0; <a <a, 0 FI}, 
-1 本 
其 中 a = (jp,0), p1 < pa ol<oa poj(z)=(V3ro) ec- /ao 。 记 
A= {(u,0): p1 SHS p2, u#0;0<01<0<02,0F 1}, 


则 集合 T(&) 就 含有 无 穷 多 个 正 态 分 布 。 

车 过 程 发 生 异 变 后 ， 我 们 能 很 快 地 诊断 出 它 所 服从 的 概率 分 布 ， 这 当然 求 之 不 得 。 但 
实际 上 ， 要 在 无 穷 多 可 能 的 概率 分 布 中 迅速 地 推断 是 哪 一 个 概率 分 布 ， 是 非常 困难 的 。 
因此 ， 我 们 退 而 求 其 次 ， 看 能 否 通过 构造 多 个 控制 图 来 推断 过 程 发 生 异 变 后 的 概率 分 布 
属于 集合 FT(&0 的 某 个 小 区 域 ( 子 集合 )。 显 然 ， 区 域 越 小 ， 越 能 准确 地 推断 过 程 发 生 异 变 
后 所 服从 的 概率 分 布 。 基 于 这 一 想法 ， 我 们 在 监测 过 程 之 前 ， 就 预先 在 T(20 PER mA 
互 不 相同 的 元 素 ，pa, 1 <k <m, Fik Kullblak-Leibler fa AHF R KTA > Mm 44} Hl 
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以 pao, Sk mm 为 “中 心 点 ”的 互 不 相交 的 小 区 域 。 这 里 ,我们 称 po,, 1 <k smm 
参考 分 布 。 两 个 概率 分 布 po 和 pw 的 Kullblak-Leibler 信息 距离 了 (pa, pao) 定义 好 下 


IT(pa, Da!) = Ea | log Pa (X) 
其 中 ， 随 机 变量 X 服从 概率 分 布 pze， 而 Ea KRDM pa 所 对 应 的 数学 期 望 。 设 
Je = {Pa € TA) : Tpas Par) < min Ipa» Pa;) }, (5) 
且 令 
T=, Tk = Je- InN (UR Ji), 


HP2<k<m X, TREL, k=1,2 ,m， 互 不 相交 并 且 URe_ Tk = r(A. HE 
义 (5) WS, Fp. € Iro WIX Kullblak-Leibler 信息 距离 po 等 于 或 靠近 po o 

对 于 取 定 的 po, lS k < m， 我 们 可 构造 由 m 个 CUSUM 控制 图 所 形成 的 m 重 CUSUM 控 
WR Tmo KEX F l 


Tmc = feces {T(ar, cr) }, 


其 中 


Par (X(t) | 
Pao(X(t)) J” 


ck > 0, 1 <k<m 是 m 个 上 控制 线 。 下 面 我 们 就 利用 这 个 多 重 CUSUM 控制 图 Tmc 对 过 程 的 
异 变 进行 诊断 。 为 此 ， 我 们 先 给 出 多 重 控制 图 的 诊断 能 力 的 定义 。 
设 焉 是 一 个 二 维 的 概率 侧 度 ， 它 可 表示 为 


P(A x A’) = f Q(da)P,(A’), AEA, AEF, 
A 


n 
T (Qk, ck) = min {n : max, 》 Ze > ck}, z® = log | 
<isn = 


其 中 4 和 大 分 别 是 集合 % 和 集合 有 的 - 域 ，@ 是 定义 在 (21, A) 上 的 一 个 概率 测度 ，(Q, 下 ) 是 
随机 过 程 的 概率 空间 ， 而 Pa 是 由 pe 所 确定 的 概率 测度 。 令 2 = {a :pe Te}, 1 <k <me 
一 个 凤 重 控制 图 Tm 的 诊断 能 力 ，DC(T)， 定 义 如 下 


DC(Tm) = min { P(pa € Lr |Tm = T(@k,ck)) }, (6) 


其 中 
Jax, Q(da)Pa(Tm = 工 (ak， ck)) 
7 Ja Q(da)Pa(Tm = T(ak, Ck)) ` 


由 此 定义 ， 当 DPC(Tm) = 1 时 ， 等 价 于 P(ps € Tk|Tm = T(wk,ck)) = 1,，1<k<m， 而 这 
意味 着 一 旦 了 ,= T(wi,ckp)， 即 TT(wn,cx) 首先 发 出 过 程 变化 的 警示 ， 则 必 有 pu 属于 Te， 也 就 
是 说 过 程 变化 后 的 概率 分 布 或 为 po 或 者 与 pu, 很 靠近 。 

在 文献 [30] 中 ， 我 们 证 明了 当 工 = EQ, (T(ak,ce)) (1 <k< m) 趋 于 无 穷 时 ， 多 重 CUSUM 
控制 图 Tc 的 诊断 能 力 趋 于 1， 即 jim DC(Tmc) = 1。 这 表明 多 重 CUSUM 控制 图 的 诊断 能 
力 可 以 渐进 达到 最 大 值 。 自 然 ， 我 们 会 进一步 问 : 


P(Pa € Tk | Tm = T (ax, ck)) 
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问题 3 多 重 CUSUM 控 制图 Znc 的 诊断 能 力 是 否 是 最 优 的 ? 即 是 否 有 DPC(Tnc) = 
maxr, {DC(Tm)}? 

根据 多 重 控制 图 诊断 能 力 的 定义 我 们 知道 ， 诊 断 能 力 不 仅 与 如 何 构造 控制 图 有 关 ， 而 且 还 
与 选取 什么 样 的 参考 分 布 ， 选 多 少 参考 分 布 以 及 参考 分 布 在 空间 下 (& 中 所 处 的 位 置 都 有 密切 
的 关系 。 一 般 来 说 ， 参 考分 布 越 多 诊断 的 精确 性 越 高 ， 但 误 判 率 也 随 之 增加 。 因 此 ， 问 题 的 关 
键 是 : 

问题 4 建立 什么 标准 来 取 仿 参考 分 布 ? 


4 变 点 时 间 估 计 问 题 


监测 一 个 随机 过 程 ， 除 了 要 尽 可 能 快 地 知晓 它 是 否 发 生 了 异常 变化 以 及 异常 变化 的 大 小 、 
程度 和 范围 ， 还 需要 及 时 地 推断 过 程 是 在 什么 时 间 发 生 异 常 变化 的 。 因 为 ， 知 道 了 过 程 发 生 异 
常 变化 的 时 间 ， 将 有 助 于 我 们 找到 过 程 发 生 异 变 的 原因 。 因 此 ， 利 用 被 监测 过 程 的 统计 数据 对 
过 程 发 生 异 常 变化 的 时 间 进 行 统 计 估 计 和 分 析 一 称 之 为 变 点 时 间 估 计 ， 也 是 SPC 的 一 个 重要 
问题 。 有 关 这 个 问题 的 最 新 研究 进展 ， 可 参阅 文献 [31-33]。 这 里 我 们 着 重 介 绍 变 点 时 间 估 计 的 
三 种 常用 的 方法 ， 参 见 文献 [34,35]。 

4.1 最 大 似 然 估计 

设 相 互 独立 的 随机 过 程 {X : t = 12…:} 在 变 点 时 间 r 之 前 服从 概率 分 布 (或 概率 密 
BE) po = p(zi ao)， 而 在 7 时 刻 之 后 变 为 分 布 pi = p(zi al)， 其 中 ao 和 a 是 两 个 参数 。 当 1 < 
Tn 一 1 时 ，(X1,… ,Xn) 的 联合 概率 分 布 (或 概率 密度 ) 为 


L(ao, 1) = [I pesia) ll P( Tk; a1). 
大 一 T 十 1 
# ao 和 ai 已 知 ， 则 变 点 时 间 了 的 最 大 似 然 估计 议 可 表示 为 
和 二 arg max [Z erua) D logp (Xi;o) |, 


l<i<n 
k=i+1 


EF arg 表示 使 得 方 括号 中 的 数值 达到 最 大 所 对 应 的 ?。 若 ao Ha, RA, S Goi Mai. 分 别 为 
在 7 二 i 的 条 件 下 关于 ao 和 al 的 最 大 似 然 估计 ， 则 变 点 时 间 7 的 最 大 似 然 估计 先 可 写 为 


Pa = arg MAX ， | Sloep (Xd) + 5 log p (Xr; âri) |- 
k=i+1 


4.2 Schwarz 信息 判别 (SIC) 估计 
R L(Go:, G1i) 表示 在 7 = i WARE REF ao 和 al 的 最 大 似 然 函 数 ， 则 Schwarz 信息 判别 统 
计量 为 


SIC(i) = —2 log L(Goi; Qi) + 2logn. 


特别 当 po = p(x; ao) Fl py = p(x; al) 分 别 是 正 态 分 布 Ni 02) W N (u, 0?) 的 密度 函数 时 ,SIC() 
可 写 为 


SIC(i) = nlog 2m + ilog &, + (n ~i)loga?;+n+2ilogn, 
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其 中 为 已 知 常数 


=1 Te PT! 


ee 去 和 
ai = D(X —p), â= 》、 (Xk — u’. 


利用 Schwarz 信 息 判 别 统计 量 SIC(i)， 变 点 时 间 7 HY Schwarz 信息 判别 估计 量 人 车 , 可 表示 为 
Tn = arg ganin_{SIC(#)}. 


4.3 Bayesian 估计 

为 叙述 简便 ， 假设 {Xi; :7 一 1,2,---} 相互 独立 ， Ai ~ N(pai, 1), FFA yy = = pr F 
Urti 三 … = ne FBT, Un AA = jr -AHr+l 为 三 个 相互 独立 的 随机 变量 ， 并 且 P(r = i) = 
pli), Un ~ N(0,07), A ~ N(0,67). Ma, 8 充分 大 时 ， 变 点 时 间 r 的 后 验 分 布 可 近似 地 表示 
为 


P(r =i] Xn , Xn) = — — 


Hpi<i<n-1, 
wor p(ije™/? T= (n —%)71 Didy Xk i! Flai Xr 


Bin- | Vi-i+(n—i)-i 

我 们 知道 ， 对 于 相互 独立 的 随机 过 程 {X， : t = 1,2,…}， 其 变 点 时 间 7+ 的 最 大 似 然 估 
计 先 具有 强 一 臻 性， 并且 萝 ,不 收敛 于 7 的 概率 随 着 nw 趋 于 无 穷 是 按 指数 衰减 的 ( 常 称 为 概率 收 
敛 速率 )。 但 是 ， 一 个 真实 地 随机 系统 其 运动 过 程 大 多 是 不 独立 的 。 那 么 ， 

问题 5 当 所 监测 的 过 程 不 独立 时 ， 化 如 它 是 Markov 过程， 其 变 点 时 间 的 最 大 似 然 估 
WT, 是 否 是 一 致 的 ?概率 收敛 速率 是 什么 ? 

上 述 我 们 都 是 在 大 样本 的 前 提 下 讨论 随机 过 程 的 监测 、 诊 断 和 变 点 时 间 的 估计 问题 。 而 在 
线 监测 时 ， 通 常 我 们 得 到 的 统计 数据 都 是 小 样本 。 因 此 ， 最 亟待 解决 的 问题 是 ， 

问题 6 对 于 小 样本 统计 数据 ， 我 们 能 否 将 过 程 的 监测 、 诊 断 和 变 点 时 间 估 计 3 个 问题 统一 
地 进行 处 理 并 得 到 最 优 的 解决 方案 ? 


5 ”结束 语 


本 文 粗浅 地 介绍 了 SPC 的 三 个 基本 问题 和 研究 方法 并 附带 提出 了 6 个 吸 待 解决 的 问题 。 
实际 上 ， 关 于 SPC 的 理论 研究 成 果 非 常 丰富 ， 未 解决 的 理论 问题 也 很 多 。 想 要 详细 地 了 
解 SPC 的 理论 研究 进展 ， 可 见 文献 [35-43]。 若 对 SPC 的 应 用 有 兴趣 也 可 参见 文献 [1,44-46]。 
真诚 地 希望 通过 本 文 的 介绍 能 有 更 多 的 人 了 解 和 认识 SPC。 
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Three Basic Problems in Statistic Process Control 
HAN Dong!, TSUNG Fu-gee? 


(1- Department of Mathematics, Shanghai Jiaotong University, Shanghai 200240; 
2- Department of Industrial Engineering and Logistics Management, Hong Kong University 
of Science and Technology, Hong Kong) 


Abstract: Based on the three basic concepts: detection, diagnosis and change-point estimation in the 
statistic process control, we first introduce the optimal control chart, the CUSUM control chart, in 
detecting the distribution change of a process, then describe how to make a statistic diagnosis by using 
the CUSUM multi-chart scheme, and finally investigate three ways of estimating the change-point: 
maximum likelihood, Schwarz information criterion and Bayesian approach. Moreover, six unsolved 
problems of detection, diagnosis and change-point estimation are proposed in the paper. 

Keywords: change point detection; statistic diagnosis; quality management; statistical process control 
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